HIP-1998-18/ph 



B-ball Baryogenesis and the Baryon to 
Dark Matter Ratio. 

Kari Enqvist^ 

Department of Physics and Helsinki Institute of Physics, 
P.O. Box 9, FIN-00014 University of Helsinki, Finland 

AND 

John McDonald^ 

Department of Physics, P.O. Box 9, FIN-00014 University of Helsinki, Finland 

February 1, 2008 



Abstract 

We demonstrate that B-ball decay in the MSSM can naturally solve the puz- 
zle of why the densities of baryons and dark matter in the Universe are similar. 
This requires that the B-balls survive thermalization and decay below the freeze- 
out temperature of the neutralino LSP, typically 1-10 GeV. It is shown that this 
can happen if the baryon asymmetry originates from a squark condensate along 
the d=6 u'^(F<F D-flat direction of the MSSM scalar potential. For this to work 
the reheating temperature after inflation must be no greater than 10^~^ GeV. 
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1 Introduction 



Observations of the dynamics of galaxies and clusters of galaxies, as well as theoretical 
models of galaxy formation from primordial density fluctuations, all point to a Universe 
mainly composed of non-baryonic dark matter [|I|]. Given that, in most models of 
baryogenesis, the physics of dark matter is completely unconnected with the physics 
of baryogenesis, it is remarkable that the ratio of matter in baryons is so close to 
that in dark matter, 0.5 to 10 % in a flat Universe When one considers that 
the dark matter particles in weakly interacting cold dark matter models, such as 
super symmetric models with unbroken R-parity 0, are often heavier that the nucleons, 
this points to an even closer coincidence when expressed in terms of number density. 
Unless one can produce a convincing anthropic principle explanation of the coincidence 
of the baryon to dark matter ratio (such as could naturally occur in axion models 0), 
this strongly suggests that the origin of baryon number and dark matter are closely 
connected. 

Electroweak baryogenesis in the Minimal Supersymmetric Standard Model (MSSM) 
and its extensions 0, ^, 0] has some possibility of a connection, in the sense that both 
baryon number and dark matter have their origin in weak interaction physics. How- 
ever, the connection between the physics of dark matter freeze-out and that of 
anomalous B+L violation during the electroweak phase transition is at best tenuous, 
allowing a wide range of possible baryon to dark matter ratios. In this paper we 
will discuss an equally plausible model for baryogenesis, requiring no new low energy 
physics beyond the MSSM, which has the great advantage of being able to account for 
the baryon to dark matter ratio naturally, namely B-Ball Baryogenesis (BBB) [||. This 
is closely related to Affleck-Dine (AD) baryogenesis |^, in which a B violating scalar 
condensate forms along a D-flat direction of the MSSM scalar potential composed of 
squark and possibly slepton fields. The difference in BBB is that the condensate is 



naturally unstable with respect to the formation of Q-balls p!0| , pJ] , |T2| of baryon num- 
ber |TB[. (In the following we will use the term B-ball to refer to any Q-ball carrying 
baryon number). The condensate therefore breaks up into a mixture of B-balls and 
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free squarks carrying baryon number. 

The subsequent evolution of the B-balls will depend crucially upon the scalar po- 
tential associated with the condensate scalar, which in turn depends upon the SUSY 



breaking mechanism. One possibility being vigourously pursued [13, n, 15, 16 1 is that 
SUSY breaking occurs at low energy scales, via gauge mediated SUSY breaking [|T^. 
As a result, for field values larger than the mass of the messenger fields, the scalar po- 
tential of the condensate scalar is completely flat, resulting in B-balls whose energy per 
unit charge is proportional to B^^^^ |T^. For large enough i?, the B-ball cannot decay 
to the lightest B-carrying fermions (the nucleons) and so is completely stable. Stable 
B-balls could have a wide range of astrophysical |jl5|, [1^, |13|, experimental [1^, 0] and 
practical implications. 

A very different scenario, which is the subject of the present paper, emerges if 
SUSY breaking occurs via the more conventional supergravity hidden sector breaking 
mechanism ||T^. In this case the potential is not flat, but nevertheless radiative cor- 
rections to the (/)^-type condensate potential allow B-balls to form ||^. However, since 
the energy per unit charge in this case is roughly independent of the B-balls are 
unstable with respect to decay to quarks or nucleons. In a previous paper |^ we noted 
that the B-balls can decay at temperatures less than that of the electroweak phase 
transition. Tew, and that the decay of such B-balls could have important implications 
for baryogenesis. For example, they could protect the B asymmetry from the effects 
of L violating interactions above Te^,, when anomalous B + L violation is in thermal 
equilibrium, or allow a baryon asymmetry to come from a. B — L conserving condensate 
One goal of this paper is to consider in more detail the physics of B-ball formation, 
thermalization and decay. 

When the B-balls decay at temperatures below Tg^, the observed baryon number 
will be a combination of the baryon number originating from the decay of the B-balls 
and that from the free squarks left over after the break up of the squark condensate. 
Since the B-balls are composed of squarks, when they decay they will naturally produce 
a number density of neutralinos of the same order of magnitude as the number density 
of baryons. Therefore, if the B-balls decay sufficiently below the freeze-out temperature 
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of the lightest supersymmetric particle (LSP) neutralinos and if the number density 
of thermal relic neutralinos is less than that from B-ball decay, then the dark matter 
density and baryon density in the Universe will be naturally related by the similarity 
of their number densities. 

The actual ratio of baryons to dark matter will essentially be determined by two 
variables: (i) the mass of the neutralino LSP and (ii) the proportion of baryon number 
trapped in B-balls to that in free squarks after the break-up of the AD condensate 
(which we will refer to as the efficiency of B-ball formation). We will attempt to show 
that the baryon to dark matter ratio expected in BBB is typically of the form observed 
in the Universe, i.e. that it is naturally less than 1 but not very much smaller than 1. 

We first give an expression for the baryon to dark matter ratio in BBB. A B-ball 
will contain a B asymmetry in the form of a squark field. When the B-ball decays, for 
each unit of B produced, corresponding to the decay of 3 squarks to quarks, there will 
be at least three units of R-parity produced, corresponding to at least 3 neutralino 
LSP's. (Depending on the nature of the cascade produced by the squark decay and 
the LSP mass, more LSP pairs could be produced). Let 3 be the number of 

LSPs produced per baryon number. Let /b be the fraction of the total B asymmetry 
which is contained in the B-balls. Then the baryon to dark matter ratio will be given 
by 

pB rrin ,^ s 

PDM N^tBm^ 

where is the nucleon mass and is the neutralino LSP mass. We first suggest that 
it is rather natural to have < 1. The present LEP lower bound on the neutralino 



mass in the MSSM, assuming no constraints on the scalar masses, is 17 GeV [|1^]. If 
we were to assume radiative electroweak symmetry breaking and universal masses for 
the squarks and Higgs scalars at the unification scale, then the lower bound would 
become > 40 GeV for tan(3 < 3 ig. For A^"^ > 3, and with > 17 (40) GeV, we 
find that < 1 occurs for fs ^ 0.02 (0.008). Thus so long as more than 2% of the 
baryon asymmetry is trapped in B-balls, the observed dominance of dark matter in 
the Universe will be naturally explained. We will show that this is likely to occur. 
Primordial nucleosynthesis 0, bounds the density of baryons in the Universe 



to satisfy 0.0048 < Qsh^ ^ 0.013, where OA < h < 1. (We adopt the bound based 
on "reasonable" hmits on primordial element abundances [Q). Thus the observed 
baryon to dark matter ratio, rs ~ — ^b) (assuming a fiat Universe), satisfies 

0.005 ^ tb Iz 0.09. This can be accounted for by BBB if 

3.7 GeV < fsm^ < 67 GeV . (2) 

For example, if the LSP mass satisfies 17 (40) GeV ^ ^ 500 GeV, then the observed 
baryon to dark matter ratio can be satisfied by a wide range of fs, 0.007 'z fB{N^/3) 
^ 3.9 (1.7). We also note that if m^^, ^ 67 GeV then we must have fs < 1, implying 
that the observed baryon asymmetry must come from a mixture of decaying B-balls 
and free baryons. 

This all assumes that the asymmetry not trapped in the B-balls can survive down 
to temperatures below Tg^. However, if we were to consider a, B — L conserving 
condensate or additional L violating interactions in thermal equilibrium above Te^;, 
then the only B asymmetry which could survive anomalous B+L violation is that 
associated with the B-balls. In this case, fs would be effectively equal to 1 (we refer 
to this case as "pure" BBB) and so would have to be less than 67 GeV. 

A crucial assumption in all this is that there is effectively no subsequent annihilation 
of the LSP's coming from B-ball decays. If this were not the case, then we would lose 
the similarity of the number densities of baryons and dark matter particles which is 
essential if we are to explain the baryon to dark matter ratio naturally via B-ball 
decays. 

In order to find out if BBB can naturally account for the baryon to dark matter 
ratio in the MSSM, we must consider in some detail the formation of B-balls from 
the primordial AD condensate, their survival in the thermal environment following 
reheating and their eventual decay to baryons and dark matter. The paper is organized 
as follows. In section 2 we discuss the efficiency of B-ball formation in a minimal 
cosmological scenario. In section 3 we consider the condition for the B-balls to survive 
thermalization down to temperatures less than that of the electroweak phase transition. 
In section 4 we discuss the decay of the B-balls and the resulting LSP dark matter 
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density. In section 5 we apply our results to realistic examples. In section 6 we give 
our conclusions. In the Appendix we discuss some aspects of B-ball solutions of the 
scalar field equations. 



2 B-ball formation 

We will assume that the B asymmetry is conserved from the time of AD conden- 
sate formation until the present. In order for the baryon to dark matter ratio to be 
explained by B-ball decay in the MSSM it is necessary that B-balls can be formed 
efficiently. This means that greater than around 2% of the B asymmetry present in 
the original condensate should be trapped in the B-balls. B-balls form because of 
the attractive force due to the logarithmic radiative correction term in the condensate 
scalar potential 

where d is the dimension of the non-renormalizible term in the superpotential. Of 
particular interest to us here will be the d=6 u'^d^d'^ "squark" direction, with a non- 
renormalizible superpotential term of the form {u'^d'^d'^Y, and the d=4 u'^u'^d'^e'^ direc- 
tion, which conserves B — L. (F- and D-flat directions in the MSSM are classified in 
reference |2^). The logarithmic term causes any space-dependent perturbation in the 
condensate field to grow and go non-linear at some time, leading to the formation of 
B-balls. 

The magnitude of K will be important in our numerical estimates. From the 1-loop 
effective potential [^], for the u'^d'^d'^ direction, the correction due to gauginos with 
SUSY breaking masses is given by 

K^-- E (4) 

a, gauginos 

where the sum is over those gauginos which gain a mass from the condensate scalar 
0. The main contribution will come from the three gluinos which gain masses from 
the squark expectation values. With ~ 0.1 we obtain \K\ ^ 0.004(M3/m)2. Thus, 



depending on the ratio of the SUSY breaking gluino mass to the squark mass, we 
expect \K\ to be typically in the range 0.01 to 0.1. 

2.1 Cosmological scenario 

In one particularly likely cosmological scenario, the AD scalar obtains a negative order 
correction to its mass squared term |2^, rn? — > {rn? — cH"^) with c ~ 1, as a 



result of non-minimal kinetic terms coupling the AD scalar to the inflaton and other 
fields in the supergravity Kahler potential. After inflation, the energy density of the 
Universe will be matter dominated by the oscillations of the inflaton fleld about the 
minimum of its potential. This will continue until the inflaton completely decays, 
leaving a radiation dominated Universe at the reheating temperature Tr, where Tr 
must be less than about 10^~^ GeV for 7713/2 in the range 100 GeV to 1 TeV in order 
to avoid thermally regenerating gravitinos 0]. In this scenario there is an initial 
spectrum of perturbations of the magnitude of AD fleld, due to quantum fluctuations 
during inflation, given by [§ 

Here Hi is the Hubble parameter during inflation and Aq is the length scale at ~ iv? . 
This is the spectrum of perturbations when ^ m^, at which time the effective mass 
squared term changes sign and the AD fleld starts oscillating, forming the condensate. 
(In the following we will use the convention that subscript o denotes the time at which 
the AD condensate starts oscillating and i denotes the time when the perturbations 
forming the B-balls go non-linear). 

2.2 Growth of perturbations 

In order to discuss the linear evolution of the perturbations, it is convenient to consider 



the special case of a homogeneous condensate of the form |13, 22 



$ = ^e^^(*) , (6) 
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where 0(t) = {ao/a)^^'^(l)o and 0{tY ~ m^, up to corrections of order Km^, and where 
we assume that \K\ is small compared with 1. This describes a condensate with the 
maximum possible charge asymmetry for a given maximum amplitude of $. Since it 
may be shown that the growth of perturbations during the linear regime is unaffected 
by the charge of the condensate, the results obtained for this homogeneous condensate 
(which is particularly easy to perturb about) will also apply to the more general case 
with a smaller charge asymmetry. 

The linearized equations for the perturbations (j) = (j){t) + 6(j){x,t) and 6 = 6{t) + 
S9{x,t) are given by [1^ 



y2 

64> + 3H6^ - {2d{t)(P{t)5d + 6<pd{tf) -5(p = -U"{^{t))64> (7) 



a2 



and 

(j){t)se + 3H{e{t)6(j) + se(f){t)) + 2{(j){t)6e + 6<i)e{t)) - ^V'5e = o 



Assuming a solution of the form |T^ 



3/2 



/,^gi(5(t)+k.x) 




a 

69 = 5^„e'(^W+'^-") (9) 

where 6(f)o and S9o are the initial values of the perturbations, and, solving for S{t), the 
perturbations grow according to 

\2a^ 9{tf ) j ^ ' 

and 

These apply if |k^/a^| ^ \2Km?\, is small compared with and l-f^l ^ 1. If 
the first condition is not satisfied, the gradient energy of the perturbations produces 
a positive pressure larger than the negative pressure due to the attractive force from 
the logarithmic term, preventing the growth of the perturbations. 
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2.3 Evolution of the condensate lumps 

The initial size and charge of the non-hnear lumps of condensate will be determined 
by the scale of the first perturbation to go non-linear. When the lumps initially go 
non-linear, they will be roughly spherically symmetric. We will model these lumps by 
a homogeneous condensate field of the form 

$ = A{t)Cos{mt) + iB{t)Sin{mt) ; r < R{t) 

$ = ; r > R{t) , (12) 

where A{t) 3> B(t). (From now on we will consider a condensate with a small charge 
asymmetry). Non-zero B{t) gives the condensate a charge, 

Q = -iJ d^x{^^^ - $t$) = ^R^mAB , (13) 

where in the final expression we have assumed that AB — AB, which will be shown 
to be consistent. 

Physically, we can understand how the lump will evolve. Roughly speaking, the 
particles in the condensate lump, which we take to be initially at rest, will be attracted 
towards each other by the force from the logarithmic term. Thus the radius of the 
condensate lump will decrease, and the value of A{t) will increase. This will continue 
until the value of A becomes so large that the non-renormalizible term in the scalar 
potential or the gradient energy dominates the logarithmic term, at which point the 
force between the particles in the lump will become rapidly repulsive and the lump will 
"bounce" elastically and start to expand. In the absence of a mechanism to dissipate 
the energy, the lump would simply pulsate. We consider first the condition for this 
pulsating lump to be able to retain a given charge. Assuming that the lump remains 
homogeneous throughout (there could be formation of a thin- wall of width of the order 
of m~^, giving the lump the profile of a thin- wall bubble, but this will not affect the 
dynamics significantly), the energy of the lump, up to terms proportional to \K\, will 
be given by E Ei + Eg, where 

E= I d^x + |V$p + U{^) , (14) 



is the total energy of the condensate lump, 

El ^ ^RVA^ (15) 

is the energy in of lump in the absence of charge and 

Eq ^ ^R^m^B^ (16) 

is the energy associated with the charge. Since the energy of the lump is constant as 
it! decreases, we have A{t) oc R"^/'^. If the charge is trapped in the lump and so is 
conserved, we will have B{t) oc 1/{AR^) cx (This implies that AB = AB, as 

we have assumed). Thus the energy associated with the charge (oc R^B) will remain 
constant as the lump contracts. Therefore if the energy per unit charge of the initial 
lump is small enough for the charge to be trapped in the lump. Eg < mQ, then it will 
remain trapped throughout the pulsation of the lump. This condition is satisfied if 

Q < Qmax = ^mR'A^ . (17) 

Thus if we consider the initial condensate to have a charge asymmetry density corre- 
sponding to the baryon asymmetry of the Universe, ub, we see that the whole charge 
within the initial lump can be trapped if initially nsi ^ AmAf, where Ai is the initial 
amplitude of the scalar field and n^j is the baryon asymmetry when the lump first 
goes non-linear. 

In general, the energy of the initial lump will be larger than that of the final B-ball 
of the same charge. Thus there must be some mechanism by which the lump can lose 
energy. For example, this could occur by the radiation of classical scalar field waves 
created by the pulsation of the lump. We will not attempt to discuss the details of 
the energy loss mechanism here, since our arguement is based purely on the fact that 
it is energetically favourable for the charge to remain within the lump. The result of 
such an energy loss mechanism is that the value of R and A when the lump reaches 
its largest size will decrease with time. However, by charge conservation, since the 
charge is proportional to ABR^, the value of B will increase. Eventually the values 
of A and B will become similar, at which point the charge of the lump will dominate 
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its subsequent evolution. The B-ball will form soon after this point is reached. The 
smallest possible values of R and A during the evolution of the condensate lump will 
correspond to the values associated with a B-ball of the corresponding charge, Rb and 
Ab- Thus so long as the condition B ^ B^axi^B, Rb) is still satisfied, we expect that 
the B-ball will be able to form from the collapsing condensate lump without any loss 
of charge. 

In general, a B-ball can either be thick or thin- walled. In the case of the thin-walled 
B-ball, the radius is a function of charge whilst the value of A remains constant 

Ab = A,; Rb = CB'/\ (18) 

where Ac is the value of the field inside the thin-wall B-ball (we discuss some properties 
of B-ball solutions in the Appendix). A thin- walled B-ball will form if the charge is 
greater than some critical charge. Be, which depends on details of the scalar potential, 
in particular on its dimension d. In the case of the thick-walled B-ball, corresponding 
to a charge less than B^, the radius remains approximately constant whilst A becomes 
charge dependent 

Ab^vB^/^; Rb^Rc. (19) 

In both cases, we see that the condition for a lump of a given charge B to be able to 
form a B-ball without losing charge, B < Bmaxi^B, Rb), reduces to a B-independent 
condition (which we will refer to as the B-ball condition) ; for the thin- walled case, 

1 < ^-fm^Al (20) 

whilst for the thick-walled case, 

1 < ^-fmRln' . (21) 

Therefore we can conclude that if the charge of the initial lump satisfies the condition 
to be trapped in the lump {jib "Z 4:mAf) and if the appropriate B-ball condition is 
satisfied, then B-balls may be expected to form with a high efficiency, with a large 
percentage of the charge of the condensate being trapped in the B-ball. 
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In fact, we can show that the B-ball conditions will be generally satisfied for any B- 
ball. We first review some properties of thin and thick- walled B-balls (see Appendix). 
In general the B-ball is described by fl^ 



(23) 



$(r,t) = ^e-S (22) 

where 0(r) and the constant uo are found by minimizing the energy of the scalar field 
configuration for a fixed charge. For the thin-walled B-ball, the energy per unit charge 
and radius are given by JTU 

E_ (2U{(t),^^^'^ 
B 

and 

R={^ ^-^V- (24) 

V4vr(202f/(0j)V2y 

From numerically solving for the B-ball for \K\ in the range 0.01 to 0.1, we find that 
the thick-wall B-ball radius is given by ~ kii{\K\^/'-^m)~^ ^ where ^ 1.2 (see 
Appendix). (From now on we will use Rc ~ {\K\^/'^m)^^). The critical charge Be 
corresponds to the charge at which the thin- wall radius, Eq. (|2^), is approximately 
equal to i?c, 

The field inside the thin- walled B-ball, corresponds to the minimum of E/B as 
a function of (j). The value of C is then set by the values of R and B at which the 
thin- wall B-ball reaches the thick- wall limit, 

/ 3 \ 1/3 X 

The thin-walled B-ball condition then becomes 

, 2mB , , 

1 < , (27) 

where we have used = = Since for the B-ball the energy per unit charge, 

E/B, is always less than m, this condition will be generally satisfied. Similarly, for the 
thick-walled B-ball, we can fix the value of rj when the charge equals Be and the field 
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at the centre of the thick-walled B-ball equals Ac- The thick-walled B-ball condition 
is then found to be the same as for the thin-walled case. Thus we can conclude that 
the B-ball condition will be generally satisfied. 



2.4 B-ball charge 

The initial charge of the B-balls will be determined by the wavelength of the first 
perturbation mode to go non-linear. This will correspond to the first perturbation 
which can result in 50/0 ~ 1 in a time of the order of H~^{t). For an inflaton 
oscillation matter dominated Universe, the exponential growth factor is 

(We take the scale factor when the AD oscillations begin to be equal to 1). The largest 

growth factor will correspond to the largest value of for which growth can occur, 
k^/a^ ~ 2\K\rf? . Thus the value of H at which the first perturbation goes non- linear 
will be 

Hi ~ -Vy , 29 

with 

a{K) = -log , (30) 

where Aq is the length scale of the perturbation ai H ^ m and 0o is the value of (p 
when the AD condensate oscillations begin. The charge of the condensate lump is 

determined by the baryon asymmetry of the Universe at Hi and the initial size of the 
perturbation when it goes non-linear. The initial non-linear region has a radius Aj at 
Hi given by 

TT 

^ {a{Xo)mH,y/^ ■ ^^^^ 

The baryon asymmetry of the Universe at a given value of H during inflaton oscillation 
domination is given by 



riB 



12 



where we have taken the baryon to entropy ratio to be rjB ~ 10 Thus the charge 
in the initial condensate lump is given by 

P 47r3 Vb\K\'/'M'^, ini5|;riV2 /100GeV\ GeV \ MO^ / vb 

"" = ^2 ma^Tn = ' ^ '° l"^' l^^^ji^^JUj llO^ 

(33) 

where we have used a(Ao) = 40 as a typical value. It is important to note that 
the details of the initial spectrum of perturbations and the condensate potential only 
enter logarithmically through a{\o)- The condition for the B-balls to form efficiently, 
Eq. ([T7|), may then be written as 



> VBmMl, _ ( m \ / 4.1xl0^^ GeV y / r^s \ . . 

where = 4.1 x 10^^ GeV is a typical initial value for the d=6 u'^d'^d'^ direction. We 
will apply these results to realistic examples later. 



3 B-ball thermalization 

We next consider the conditions under which the B-balls can survive down to temper- 
atures less than that of the electroweak phase transition Tf,yj. In order to discuss this 
question, we will consider throughout a Gaussian Ansatz for the thick-walled B-ball 

0(r) = 0(O)e"^ , (35) 

where 0(0) is the value of at the center of the B-ball. In the Appendix it is shown 
that the Gaussian is a physically reasonable approximation for the thick-walled B-ball. 

There are two processes which might lead to the destruction of B-balls by the 
thermal background. The first process is the dissociation of the B-ball by collisions 
of thermal particles with the "hard core" of the B-ball, corresponding the region of 
the B-ball that thermal particles cannot penetrate. (This was considered first in 
reference |l^). The second process is the dissolution of the B-ball by the thermalization 



and transport of charge from the "soft edge" of the B-ball, within which a thermal 
equilibrium distribution of background particles can exist p[ . 
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3.1 Dissociation 

We first consider the dissociation of the B-ball. Thermal particles can penetrate the 
B-ball to the radius at which g'0(r) ST; we refer to this as the stopping radius r^t, 
with the corresponding field being (j)st- (This defines the hard core of the B-ball). As 
the thermal particles enter the B-ball, they will transfer energy to the classical scalar 
field, finally coming to a (temporary) halt. If the time scale over which the particle 
stops is short compared with the time over which the scalar field of the B-ball can 
absorb the energy, 6te ^ l/m (set by the dynamical scale of the scalar field) , then the 
thermal particle will be ejected, leaving some energy in the hard part of the B-ball, 
with the amount of energy transferred depending on how quickly the particle comes 
to a halt. Let 6tr = Kj./m {K^ ^ 1) be the time scale over which an excited B-ball can 
radiate its excess energy in the form of scalar field waves. If sufficient energy can be 
delivered to the B-ball within 5tr to overcome the binding energy of the charges in the 
B-ball, then the B-ball will dissociate. If, on the other hand, less energy is delivered 
in this time, then the B-ball will be able to radiate the excess energy adiabatically 
and so will not dissociate. Thus the condition for the B-ball to evade dissociation is 
that AE{Str) < 5mB, where l^E{5tr) is the energy delivered to the B-ball by thermal 
particles in a time 8tr and 5m — {m — E/B) ^ \K\m (see Appendix) is the binding 
energy per unit charge. 

The flux of particles onto the hard core of the B-ball is 



where g{T) fti 100 is the effective number of light thermal degrees of freedom coupling 
to the condensate scalars, and, for the Gaussian thick wall B-ball, 



The energy transferred in a collision will depend on how quickly the field of the 
B-ball can respond to the impacting thermal particle. The time scale over which a 
particle of energy 3T comes to a halt in the hard core is roughly the distance over 
which the change in the effective mass of the incoming particle, g5(f), is approximately 




(36) 




(37) 
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equal to 3T around Vst- For the Gaussian wall this corresponds to 

5rst ^ — = 4 • (38) 

Thus we expect the suppression factor in the transfer of energy to be of the order of 
m~^/5rst = 2(3 /{niR) ^ 2l3\K\^/'^^ which is typically less than 1 but not very much 
so. Therefore we do not expect a very large suppression in the energy transferred. 
(In fact, we will show later that a thermal equilibrium exists within the soft edge of 
the B-ball. Therefore thermal scattering may slow an incoming particle, so that the 
B-ball will have time to absorb most of its energy). Let the energy per thermal particle 
transferred to the B-ball be 7t7', where 7t ^ 3. Then the rate of energy increase due 
to incoming thermal particles will be 

dE 4g{T)^TT^f3^R^ 



dt 71 



(39) 



The requirement to avoid dissociation is that, in the time 6tr, insufficient energy is 
supplied to the B-ball to overcome the binding energy. This gives an upper bound on 
the temperature at which B-balls can exist: 

.1/4 



T 



n\K\ ( 5m 



mB^'^ . (40) 



_4^(T)i^,7T/32 V m 

Combining this with the expression for the B-ball charge as a function of the reheating 
temperature, Eq. (|33D , then gives an upper bound on the reheating temperature 

1 / 5m\ 



Tr < 1 X lO^lirp/io ' 



100 GeV 



J 



a) VlO-io/ 



(41) 

where we have taken g{T) ^ 100. Thus, using as typical values ~ 10, a ~ 40, 
/? 4, 5m/m ^ \K\ and 7r ~ 1, we obtain for \K\ 0.01 (0.1) 

/ <rj) \ 3/5 

Therefore, so long as Tr ^ 10^~^ GeV, the B-balls should be able to evade dissociation 
for li^'l in the range 0.01 to 0.1. 
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3.2 Dissolution 

The second process which might lead to the thermahzation of the B-ball is dissolution 
by the removal of charge from the soft edge of the B-ball. The width of the soft edge 
will correspond to the distance at Vst over which does not change much, 50/0 ^ 1. 
For the Gaussian B-ball this corresponds to 5rst ~ R/{2P) ^ m/(2\K\^^'^P). This is 
much larger than the mean free path of the strongly interacting thermal background 
quarks, Xmfp ~ kg/T (where fcg ~ 6 |^), ifT ^ m. So we expect a thermal equilibrium 
to exist within the soft edge at high temperatures. We also expect the B-ball charges 
within the soft edge to become thermalized. The rate at which charge is removed 
from the B-ball will then depend on two factors. Firstly, how rapidly the thermalized 
charges diffuse out of the B-ball. Secondly, how rapidly the hard core of the B-ball 
reconfigures itself in order to compensate for the loss of charge and so minimize its 
energy, replenishing the charge in the soft edge. The rate of diffusion may be estimated 
as follows. The number of steps for a charge to leave the soft edge by a random walk 
will be of the order of {Srst/ XmfpY ~ {RT / {2kqj3)y . The time for each step is around 
Xmfp- Therefore the time taken for a charge to leave the soft edge is 

The total charge in the soft edge is Bsoft ~ ^'^^'t'1t'^1t^''^st- Thus the rate at which 
charge can be removed, assuming that the B-ball can reconfigure itself on a time scale 
shorter that r^, is 

IdB _l Bsoft ^^^^ 



B dt B Td 

where B ^ 2uj(f){0)'^R^ is the total charge of the Gaussian B-ball. The time for the 
B-ball to reconfigure itself will be of the order of m~^, which will be short compared 
with Td for T ^ = 4j3'^kq\K\m. So at high temperatures we expect that the B-ball 
can efficiently replenish the charge lost by thermahzation of the edge. The rate at 
which charge is lost is then 

1 dB _ Airp^kgT 

B~dt ^ ^20(o)2i?2 • y ) 
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Requiring that this is less than H, in order to avoid dissolution of the B-ball, then 
imposes a lower limit on the B-ball charge 



where H = kxT'^/Mpi [kr ~ 17). For T ^ m=K, the rate at which charge is lost will 
be determined by the rate at which the B-ball can replenish the lost charge, in which 



case T in Eq. (|6) should be replaced by m* and the lower bound on B becomes T 
independent. Requiring that Eq. (^Bf) is satisfied for all T ^ then gives the lower 
bound on B, 



'P\ /lOO GeV\ 



Using Eq. (|33D we then obtain an upper bound on the reheating temperature 



Therefore, with g ^ 1, kg ^ 6 and li^l = 0.01 (0.1), we find that the reheating 
temperature cannot be larger than 10'^ (10^) GeV if the B-balls are to have sufficiently 
large charge to evade dissolution. We note that this bound is tighter than that from 
dissociation for all \K\ ^ 0.01. 

Thus we can conclude that dissociation of B-balls by thermal particle collisions 
and dissolution by soft edge thermalization can both be avoided if the reheating tem- 
perature is sufficiently low, less that about 10^ — 10^ GeV for \K\ in the range 0.01 to 
0.1. 

3.3 Dissociation prior to reheating 

The above discussion applies to the radiation dominated period following reheating. 
The "reheating" temperature in fact refers to the temperature at which the radia- 
tion from the decay of the infiaton condensate comes to dominate the matter density 
remaining in the condensate. The temperature decreases steadily during the infia- 
ton oscillation dominated period, with no "reheating" as such. The radiation energy 
density due to infiaton decays is given by fl], ^ 

P.- , (49) 
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where T is the decay rate of the inflaton condensate particles, which have a mass density 
p. Therefore we must ensure that this radiation energy density does not dissociate the 
B-balls prior to reheating. (Dissolution by thermalization of the soft edge of the B- 



ball is most effective at lower temperatures, since the rate of dissolution, Eq. (|45|), 
is proportional to T, whereas H is proportional to for radiation domination and 

for inflaton domination. Therefore the possibility of higher temperatures existing 
prior to reheating will not alter the upper bound on Tr following from dissolution). In 
addition to the inflaton condensate, there is a second source of radiation, namely the 
remains of the original squark condensate. After the AD condensate breaks up into 
B-balls and non-relativistic squarks (the initial momentum of the squarks radiating 
from the collapsing condensate lumps will be at most of the order of the inverse of the 
radius of the B-ball, \K\^^'^m), the non-relativistic squarks will still evolve as a matter 
density much like a condensate. The average value of the scalar field associated with 
this density of squarks will be of the same order as the original AD field amplitude 
would have been had the condensate not collapsed. We must also check that the 
radiation from the decays of these squarks does not dissociate the B-balls. 

We first consider inflaton decays and ignore the radiation energy coming from 
squark decays. Inflaton decays with a constant decay rate F will produce a background 
radiation energy density with a temperature given by 
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Tr ^ K{MpiHTlf'' ■ K = [-^^^^^) ' (50) 
where kj. ^ 0.4 for g{T) ^ 100. At the time when the B-balls form, corresponding to 



H = Hi, where 



H,^,.m-^(m(,^]('-^)GeV. (51) 



^0.01 ; VIOO GeV/ V a 
the temperature of the background radiation is 



The upper bound on T to avoid dissociation, Eq. (|40|) , may be written as, using 
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Eq. Q, 

2 X 10= m]'" (l^Y' C^y" GeV . (53) 

\0.01J VlOOGeV/ \ Tn J \aj ^ ' 

Thus the temperature of the background radiation is essentially low enough to evade 
the dissociation of the B-balls if Tr ^ 10^ GeV. In fact, the value of H at which the 
B-balls can be dissociated will be much less than Hi. This is because at H ^ Hi the 
squark field expectation value, due to the non-relativistic squarks from the AD 

condensate will be large enough to give an effective mass g{\4>\), much larger than T^, to 
any particles coupling directly to the B-ball squarks, so Boltzmann suppressing them 
and preventing them from dissociating the B-balls. The value of {\(f)\) as a function 
of H is given by = {H/Ho)(f)o, with H following from Eq. (^) and H^ ~ m. 
g{\4)\) ^ T is satisfied down to 

^ . 9.6 X 10- i^Y" fii2ii!?!l^V'' GeV . (54) 

VlOOGeV/ V 9(Po J VlOOGeV/ ' ^ ' 

which is much smaller than Hi for both the d=4 and d=6 condensates. (This assumes 
that the squark condensate has not decayed or thermalized; we will discuss this point 
below). The corresponding temperature is given by 

/ rn \ 1/3 /4.1 X 10^4 GeVX'^V Tn N^/^ 

Thus dissociation will not occur until H -C H^. As a result, dissociation of the B-balls 
during the inflaton dominated period by radiation from inflaton decays will impose no 
further constraint if the reheating temperature satisfies the upper bound coming from 
dissolution, Tn < lO^-^ GeV. 

In addition, we should check that the radiation from the decay of the low momen- 
tum squarks presents no problems. To show this, we neglect the radiation coming 
from the inflaton decays. The radiation from the squark decays will have a density 
given by Eq. (|^), where, for g{\(t)\) ^ T, the squark decay rate is approximately given 

by, 

r ~ yrr^ , (56) 
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where ag = 5'^/47r and where g is typically the strong gauge coupling. This decay rate 

is time dependent. As a result, the radiation temperature increases as H decreases, 

Tr oc p,y^ oc The maximum temperature of the radiation from the squark 

decays will therefore correspond to the smallest possible value of H, at which the 

energy density in radiation equals that remaining in the low momentum squarks, 

2ag\i/3 ^5/3 



^2/3 • (57) 



This gives for the maximum temperature of the radiation from squark decays 

This is generally much less than the upper bound on T coming from avoiding dissoci- 
ation, Eq. (|55|), for acceptable values of Tr. 

Finally, let us estimate the time at which the squark condensate is thermalized. 
So long as this occurs well after B-ball formation, the B-ball should be protected from 
thermalization by dissociation. Assuming that g{\4>\) ^ T, the largest thermalization 
rate possible will be of the order of 

^"i^- ^^^^ 

Thermalization will occur once T ^ H, which, assuming inflaton domination, occurs 
once H is smaller than Hthermi where 

8/9 

J 



'-therm 



(60) 



For the d=6 condensate this is much smaller than the value at which the B-balls form, 
iiTj, when Tr ^ 10^~^ GeV and \K\ is in the range 0.01 to 0.1. For the d=4 conden- 
sate, on the other hand, for which typically 0o ~ 10^° GeV, it can be larger than 
depending on the values of Tr and \K\. Thus it is possible that the d=4 squark 
condensate could be thermalized before the B-balls form. Let us finally note that the 
squark condensate will thermalize before the reheating temperature is reached (cor- 
responding to EiT^ ^ 2 X 10-i2(^Tfl/103 GeV)2 GeV) for all values of Tr consistent 
with B-balls surviving thermalization. 
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To conclude, we find that B-balls can survive thermalization both during and after 
the infiaton dominated era so long as the reheating temperature is sufficiently small, 
Tr 'Z 10^~^ GeV for \K\ in the range 0.01 to 0.1. The strongest constraint comes from 
dissolution of the B-ball by thermalization and diffusion of charge from the soft edge. 



4 B-ball decay 

Assuming that the reheating temperature after inflation is sufficiently low that the B- 
balls survive thermalization, the next question to be considered is that of their decay 
temperature. Typically the LSP freeze-out temperature is Tjr ~ mj(,/20. Therefore if 
the LSP is not very heavy, say ^ 200 GeV, then we must require that the B-balls 
decay at temperatures less than around 10 GeV, in order that they decay after the 
freeze-out of the LSPs, with the LSPs from B-ball decay dominating the relic density. 
We must also ensure that there is no subsequent annihilation of the LSPs coming from 
B-ball decays, in order to maintain the direct relationship between the baryon and 
dark matter particle number densities. 



4.1 Decay temperature 

In the previous discussions of B-ball decay, the decay rate has been estimated from 



the model of reference [0, which assumes that there is a single scalar coupled to a 
light fermion. In this case the decay to light fermions is proportional to the area of 
the B-ball (decay within the volume of the B-ball being blocked by the Pauli principle 
once the Dirac sea is filled), with an upper bound on the decay rate being given by 



dB 



dt 



fermion 

where A is the area of the B-ball. The upper bound is likely to be saturated for B- 
balls with 0(0) much larger than m. However, in realistic SUSY models, there will be 
several scalar fields and in general we would not expect that the condensate scalar will 
correspond to the lightest scalar. Since the decay to scalar fields within the volume is 
not blocked by the Pauli principle, it is possible that the decay to light scalars could 
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be significantly enhanced relative to the decay to fermions. However, the decay to 
light scalars will only be possible near the edge of the thick-walled B-ball. This is 
because particles coupling directly to the condensate scalars will gain a large effective 
mass from (0) inside the B-ball. As a result, decay to light scalars will occur only via 
loop diagrams with rates suppressed by this large effective mass. 

In the previous subsection we derived a bound on the reheating temperature from 
B-ball stability. Let us now derive a general upper bound on the reheating temperature 
from the requirement that B-balls decay below the LSP freeze-out temperature. Let 
fs be the possible enhancement factor of the scalar decay rate over the fermion decay 
rate, 

dB _ fdB\ 



fs\ — ] ■ (62) 

fermion 



dt \dt ) 

The B-balls will decay once the decay rate is larger than H . This gives for the decay 
temperature Td 

^ . I!..^R^M„ V'\, ^^^^ n\ . v. / lo^^ ^ ^^^^ 



iSnkrE / \\K\ VlOO GeV/ \B 



Using Eq. (|33|), we find that requiring that is less than 10 GeV imposes an upper 
limit on the reheating temperature 

Thus, with \K\ ^ (0.01 - 0.1), we need < (5 x 10^ - 2 x 10^)/-^ GeV for the 
B-balls to decay at a temperature below 10 GeV. 

We next consider the possible enhancement of the decay rate to light scalars over 
the decay rate to light fermions. We first note that, for the particular case of the 
d=6 u'^d'^d'^ direction, it is quite possible that there will be no scalar decay mode en- 
hancement. The condensate scalar in this case is a linear combination of right-handed 
squarks. For models with universal soft SUSY breaking terms at large renormalization 
mass scales, the solution of the renormalization group equations implies that the left- 
handed squark masses will typically be heavy compared with the right-handed squark 
masses at low mass scales [113] . The Higgs scalar masses can also be heavy compared 
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with the right-handed squark masses, depending on the jj, parameter. On the other 
hand, in models with universal soft SUSY breaking terms at the large mass scale, 
the sleptons will typically be lighter than the right-handed squarks. However, even if 
slepton masses or Higgs masses were less than the right-handed squark masses, any 
condensate squark decay to sleptons or Higgs bosons would also have to involve a quark 
in the final state, together with a gaugino-Higgsino or lepton. So this would effectively 
involve pair producing light fermions in the final state, which would be suppressed 
by the Pauli principle as for the two fcrmion decay process. The only exception to 
this would be if the mostly right-handed stop was sufficiently light compared with the 
other right-handed squarks that decay of the condensate squark to a stop plus a Higgs 
boson became kinematically possible. This is a model-dependent possibihty. There- 
fore it is quite possible that decay of the linear combination of right-handed squarks, 
which constitutes the condensate scalar in the d=6 u'^d'^d'^ direction, to pairs of lighter 
scalars will be kinematically forbidden and that fg will be effectively equal to 1. In 
this case it is sufficient that Tr be less than about 10^~^ GeV for the B-balls to decay 
below 10 GeV. 

More generally (for example, for the case of a B-ball made of left-handed squarks 
or sleptons, such as for the d=6 directions based on d'^QL or e'^LL), we would expect 
a decay mode to pairs of light scalars to exist. In this case we can estimate the largest 
possible enhancement factor using the Gaussian thick-wall Ansatz. Within the B-ball, 
for values of (p much larger than m, the lowest possible dimension operator which could 
allow the condensate scalars to decay at one- loop to light particles is the d=5 operator 



where x ^^(^ V represent the light particles and M ^ g(f), where g is the coupling of 
the heavy particles to 4>. The rate of decay of the condensate scalars to light scalars 
will then be 




(65) 




(66) 
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where a;0^(r) is the charge density within the B-ball and where 

2 'H 

am 

r(r) ^ — — ; g(j)> m 



~ am ; gcj) < m , (67) 

with a = (7^/(47r). (For simphcity we consider a single couphng constant g). Let r^, 
be the radius at which (f){r) = m/g. Then the largest contribution to the decay rate 
will come from a region of width 6r ^ i?^/ (4r*) around r*, over which has a roughly 
constant value (p ^ m/g, where 

= 7i? ; 7 = ( ^-M] . (gg) 
\^ m y 

From Eq. ( |66D this gives a rate, 

-Anaum [m'^a j drr"^ + 0^(0) j drr'^e~~R^ \ (69) 



dB 



l + ^ll^l ■ (70) 

fermion 



^_m7 / 7y\ (dB\ 
^ \K\y^ 37T )[dt) 

where we have used w ~ m (see Appendix). For the thick- walled B-ball, for typical 
values of the parameters, g(l){0)/m ~ (0.1 — 0.01)i?^/^ and so 7 ~ 4.5. Thus we find 
that the enhancement factor is typically given by 

/^-^(l + 2V) . (71) 

Thus, if g is less than 1, then for \K\ ^ 0.01 — 0.1 we expect an enhancement factor 
not much larger than about 10^. (For g less than 1 most of the enhancement factor 
comes from unsuppressed tree level decays occuring at r > r^). This would impose an 
upper bound on the reheating temperature of about 10'^~^ GeV in order to have B-ball 
decay at a temperature less than 10 GeV. Therefore, if it is possible for the condensate 
scalars to decay to light scalars, then the upper bound on Tr will be much tighter, 
although it will not be very much different from that coming from the requirement 
that the B-balls survive thermalization by dissolution. 

Thus, if Tr < 10^ (10^) GeV for \K\ = 0.01 (0.1) , then we expect that the B-balls 
will typically survive thermalization and will decay at a temperature below 10 GeV. 
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In this case we may be able to understand the baryon to dark matter ratio as being 
a consequence of B-ball decays. Just how much dark matter is produced will depend 
on the details of the thermalization and annihilation of the LSPs coming from the 
decay of the B-balls. It should be emphasized that the requirement that B balls evade 
thermalization and that they decay at a low enough temperature to account for the 
baryon to dark matter ratio both point to a relatively low reheating temperature after 
inflation. We will see that this is also preferred by the naturalness of the observed 
baryon asymmetry in this scenario. 



4.2 LSP dark matter from B-ball decay 

A typical mass of a B-ball would be 10^^ {10^ GeV/T^i) GeV. When they decay at 
t ^ the average distance between two B-balls is /b ~ 5 x 10 ^(Td/lOGeV)lH, 
where In is the horizon distance and is the decay temperature (and assuming that 
the baryon number of the Universe is mainly due to B-ball decay). As we discussed in 
the Introduction, there will be A^^^ ^ 3 LSPs produced per baryon number, so that as a 
reference number we take the total number of LSPs produced in each B-ball decay to 
be A^LSP ~ 10"^^ {10^ GeV/TR). Because the B-ball decay gives rise to an overabundance 
of LSPs with respect to the relic density, we must check that the produced LSPs will 
not get annihilated, both during the initial decay of the B-balls and subsequently once 
the LSP density is smoothed out. 

Let us consider the decay of a single B-ball. The produced LSPs will collide with 



the weakly interacting particles in the background and after a few collisions 0] will 
locally settle into a kinetic equilibrium (some radial bulk motion might still exist). 
Thus they soon become non-relativistic as is by assumption less than the freeze- 
out temperature of the LSPs, which we may take to be about 0. The rough 
freeze-out condition for LSPs initially in thermal equilibrium is 

nLSp{cr.nnv) ^ Hf-^ , (72) 

where o"ann is the LSP annihilation cross-section and the subscript / refers to the 
freeze-out values. LSPs annihilate to light fermions. The thermally averaged cross 
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section can be written as (aann'w) = a + hT /m^, where a and h depend on the couphngs 
and the masses of the hght fermions P] . In the following we will consider the case of a 
light neutralino, < mw-, and we will neglect the final state fermion masses. (This 
would be appropriate for the case of pure BBB (/^ = 1), for which ^ QTGeV). 
In this case a = and (cr^^^v) oc T. In that case one obtains from Eq. (|72D that 
b ~ Hm'^T^'^nJ^ , where 

% = -(^(^xTff'e-"^>^^^f ^ 1.46 X 10-''ml , (73) 

is the freeze-out LSP density and where we have used Tf ^ m^/20. 

The LSPs produced in the decay of the B-ball will spread out by a random walk 
with a rate u determined by the collision frequency divided by the thermal velocity 



Vth ~ yT/m^, which reads 

t^'^ ~ i{crscVrei)nrad/vth) ~ giT)G%m^T^ , (74) 

where ~ SGGjpT^/n is the scattering cross-section with thermal particles, rirad 
is the density of particles in radiation and we have taken Vj-d ~ 1. It is very likely 
that the decay is spherically symmetric. Thus the number of LSPs created at t = 
subsequently form a Gaussian distribution around the decaying B-ball with 



rfiVLSp(r,t = 0) / 2 _,./(2,t) 



dr V ^t'^ 

Folding in the rate by which the LSPs are produced, 

TB^ipe-^^' , (76) 



dt 

and integrating over time we find that the radial distrubution of the LSPs at time t is 
given by 

dN^p{r,t) _ fi I 2 tot -TBt' -ry(2u(t-t')) 

The largest contribution to the integral comes from the region t' ~ F^^ ^ u so that 
in the leading approximation Eq. (|77|) reads 



dr 



TTUX 
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where we have written t = T^^x with x ~ Thus there is a "central region" of 

radius 

'<^r (79) 

within which the mean LSP number density is given by 

'^...(.)=.5^(^1", (80) 



TTUX 



and outside of which it is exponentially suppressed. 

In the central region annihilation is significant if ^lsp (cann^rci) ^ H. Neglecting 
the possible change in the light degrees of freedom below Ty so that H/Hf = T'^/Tj, 
one finds that significant annihilation means simply that ^lsp(^) ~ /m^nf. From 
Eq. (|80|) this happens when (we now take x ~ 1) 



81) 



A large fraction of the LSPs produced in the B-ball decay will thus be destroyed if the 
critical distance is greater than the size of the central region, r = = {u/VbY^'^- 

Otherwise the build-up of LSP density is slow so that the main bulk of the LSPs will 
have time enough to escape from the vicinity of the B-ball. We find that 



r, 



GeV^ 

Thus annihilation is insignificant provided 



mv \ 3/16 / 1020 \V8 / 100 



Because it is likely that there remains some radial bulk motion even after LSP equili- 
bration, the actual limit on is bound to be somewhat less stringent. In any case, 
we may conclude that typically most of the LSPs will survive if the B-ball decay 
temperature is less than a few GeVs. 

The above discussion applies to the case of a single B-ball. However, we must also 
ensure that there is no subsequent LSP annihilation once the LSP density is smoothed 
out. This can be significant because there is typically more than one B-ball within 
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the central region of a given B-ball, so it is not sufficient to check that the LSP's can 
escape from a single B-ball without significantly annihilating. Since the essence of 
the natural explaination of the baryon to dark matter ratio via B-ball decay is the 
similarity of the number of LSPs and baryons produced by B-ball decay, we must 
ensure that the number density of LSP's is not significantly reduced by annihilations, 
as well as ensuring that the thermal relic density of LSP's is small compared with that 
from B-ball decay. 

The upper limit on the number of LSP's at a given temperature is 

niimitiT) ^ • (84) 

For an annihilation cross-section dominated by b at T^, 

niimuiT) ^ {g{Tf)/g{T)fl\Tf/Tfnr,uc{T), 
where nreiic{T) is the thermal relic density at T 

„„,„(r) « munv^) . (85) 

\9{Tt)l \Ttl \{ovUnlr, 

Thus in order to have a natural explaination of the baryon to dark matter ratio we 
must have 

nreliciTd) < riLSpiTd) < [^^^ (^) nreliciTd) , (86) 

where 

. ^ 1,3 X 10- (^)' (^) (f ) (^) GeV-= (8, 

is the LSP density from B-ball decay (assuming that the B asymmetry comes mostly 
from B-ball decay). For a reasonable range of LSP densities to exist, the B-balls must 
decay sufficiently well below the LSP freeze-out temperature. 

For the h dominated annihilation cross-section, with m^/Tf ^ 20, we find that the 
condition riLSpiTd) > nreUc(Td) is Td independent and generally satisfied, with 
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The condition that nLSp{Td) < numuiTd) is satisfied if 

-Hff" im" (^) ^-^ ■ 

Thus we can conclude, for the case of a 6 dominated annihilation cross-section, ap- 
propariate for a light neutralino, that the baryon to dark matter ratio can be naturally 
explained by B-ball decay if the B-balls decay at a temperature less than about 1 GeV. 
This in turn requires that Tr < (5 x 10^ - 2 x 10^)/^^ GeV for \K\ ^ (0.01 - 0.1). We 
note that this may favour the u'^d'^d^ direction, for which fg can be equal to 1, over 
directions along which the B-balls can decay to fight scalars, for which fg ~ 10^, so 
imposing a much tighter upper bound on the reheating temperature. 

5 Application to the u^dFdF and u^u^dFe^ directions 

In this section we apply the results of the previous sections to the u^d'^d'^ and u^u'^d'^e'^ 
squark directions. We first note that the upper bounds on the reheating temperature 
following from the requirement that the B-balls are not thermalized by dissolution and 
that they decay at a sufficiently low temperature are only very weakly dependent on 
the details of the scalar potential i.e. logarithmically via a{\o). Thus these bounds 
will be effectively independent of d. 

The initial value of the field when the condensate oscillations begin, for the d=6 
case, is given by, 

0„=5.8xlO"A-V^(^)"*GeV (90) 

and for the d=4 case by 

Thus the lower limit on the reheating temperature, following from the requirement of 
maximum efficiency for the d=6 case, is given by 
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where the natural value of A for d=6, assuming that the strength of the non-renormalizible 
interactions is set by Mp, is around 0.003. This is easily satisfied for reheating tem- 
peratures less than the upper bounds from B-ball thermalization and decay. For the 
d=4 case, the upper bound for maximum efficiency is given by 



where A is naturally around 0.1 for d=4. This cannot be satisfied if the reheating 
temperature is low enough for the B-balls to survive thermalization, Tr ^ 10^"^ GeV. 

Thus efficient B-ball formation is quite natural for the d=6 u'^d'^d'^ direction. As 
noted in the previous section, this direction is particularly favoured as = 1 is pos- 
sible, allowing a relatively large reheating temperature to be compatible with similar 
baryon and dark matter particle number densities. Therefore, so long as the reheating 
temperature after inflation is no greater than around 10^~^ GeV, the baryon to dark 
matter ratio can be naturally accounted for via B-ball decays in this case. 

For the d=4 direction, the B-balls do not form efficiently However, for the d=4 
u'^u'^d'^e'^ direction B-L is conserved. Therefore, in this case only the baryon number 
trapped in the B balls will survive. Thus it may appear that the observed B asym- 
metry can be obtained by simply having a much larger B asymmetry in the initial 
condensate than the present B asymmetry. However, for the d=4 condensate, no B- 
balls will survive thermalization if we can generate the observed baryon asymmetry. 
The maximum possible asymmetry will simply correspond to the number density of 
(f) scalars in the condensate. When the condensate scalar starts oscillating at if ~ m, 
this is given by ri^ o ~ m0^/2. The baryon to entropy ratio is related to o by 



Thus, with (/)o for d=4 given by Eq. (0), we find that the maximal baryon to entropy 
ratio initially is r/B ~ 6x 10^^°A^^(Tfj/10^ GeV) ^ 10^^. Therefore the d=4 condensate 
can only account for the observed B asymmetry if the asymmetry is not very much 
below the maximum possible for a d=4 condensate, which requires a high reheating 
temperature, Tr ^ 10^ GeV. In this case no B-balls will survive thermalization. Thus, 



Tr > 8 X lO^A GeV , 



(93) 




(94) 
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in general, the d=4 condensate will not be able account for the baryon to dark matter 
ratio via B-ball decay. 

For the case of the d=6 condensate it is important to note that, as well as allowing 
for the baryon to dark matter ratio to be explained via B-ball decay, a low reheating 
temperature is also the most natural from the point of view of the observed value of 
the B asymmetry. The observed B asymmetry for the d=6 case, assuming that B is 
conserved from the time of condensate formation until the present, is given by 



where 6c p is the CP violating phase (the phase of the A-term) responsible for the B 
asymmetry in the initial condensate. Thus we see that if Tr is small compared with 
10^ GeV then we do not need to have a CP violating phase that appears unnaturally 
small compared with 1 in order to account for the observed baryon asymmetry, rjB = 
(3 — 8) X 10^^^ 0. For example, Tr ^ 10'^ GeV would allow the present B asymmetry 
to be explained by a phase greater than or of the order of 10^^. We also note that this 
is most natural if the CP violation responsible for the B asymmetry enters through 
the A-terms in the scalar potential, rather than being set by an initial random phase 
in the AD field due to de Sitter fluctuations during inflation, which would tend to give 
Sqp ^ 1. This in turn requires that order H corrections to the A-terms exist, in order 
to damp out any large random phase, which may serve as an important constraint 
on inflation models. Thus both the baryon to dark matter ratio and the observed B 
asymmetry can be most naturally explained by having the B asymmetry originate from 
a squark condensate along a d=6 D-fiat direction of the scalar potential with a low 
reheating temperature after inflation. The u'^d'^d^ direction is particularly favoured, 
typically allowing a larger reheating temperature than the other d=6 directions in the 
MSSM by evading the large scalar mode enhancement of the B-ball decay rate. 

This all assumes that the B asymmetry from the d=6 condensate is conserved. 
In fact, this is likely to be effectively true in general, even if there were additional 
L-violating interactions, since the B-balls in this case are likely to form with a high 
efficiency, protecting most of the asymmetry. Only if the reheating temperature were 
very low, less than around 1 GeV, would there be a reduction in the efficiency of 




(95) 
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B-ball formation along the d=6 direction. However, from Eq. (0), the reheating 
temperature cannot be much less than 1 GeV and still be able to account for the 
observed asymmetry, even if all the asymmetry could survive. 

6 Conclusions 

We have considered the possibility of accounting for the baryon to dark matter ratio 
via B-ball decay in the MSSM. We have considered the constraints on the reheating 
temperature after inflation, in the context of a simple cosmological scenario, following 
from the requirements that the B-balls can survive thermalization, decay sufficiently 
long after LSP freeze-out and naturally account for the observed B asymmetry. 

For the case of a d=4 AD condensate, we find that, in general, the reheating 
temperature after infiation must be too high for the B-balls to have survived thermal- 
ization. In contrast, for the case of the d=6 u'^d'^d'^ condensate, a perfectly consistent 
cosmological scenario emerges. The requirements that the B-balls survive thermal- 
ization, can form efficiently and so naturally account for the baryon to dark matter 
ratio and can decay sufficiently below the LSP freeze out temperature without sub- 
sequent annihilation are all satisfied for typical values of the radiative correction to 
the scalar potential if the reheating temperature is between 1 GeV and 10'^"^ GeV. A 
low reheating temperature is also preferred if the B asymmetry is to originate from a 
d=6 AD condensate without requiring an unnaturally small CP violating phase. The 
u'^d'^d'^ direction is particularly favoured, allowing a larger reheating temperature than 
the other d=6 directions in the MSSM by avoiding the scalar mode enhancement of 
the B-ball decay rate. For the case of pure B-ball Baryogenesis, where all of the baryon 
asymmetry comes from B-ball decay, the LSP cannot have a mass greater than 67GeV. 

It seems remarkable that with only the fields of MSSM, together with plausible non- 
renormalizible corrections and period of primordial infiation, we are able to naturally 
account for a baryon to dark matter ratio of the observed form. Given that following 
infiation there will typically be an ensemble of domains, much larger than the observed 
Universe, in which the MSSM scalars take on all possible initial values along F- and 
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D-flat directions, it seems certain that there will be at least some domains of the 
Universe in which the scalar fields will have initial values along the u'^d'^d'^ direction. 
All we then require, in order to account for the baryon to dark matter ratio, is that 
the reheating temperature be sufficiently low. In a sense, this may be considered a 
prediction of this scenario, providing a strong motivation for infiation models with low 
reheating temperatures. The requirement of a low reheating temperature will impose a 
significant constraint on infiation models, tending to favour infiaton candidates which 
are light and/or highly decoupled from the MSSM or other light fields. 

There are several issues which remain to be discussed. The details of the dark 
matter density coming from B-ball decay and the resulting constraints on the MSSM 
parameter space should be investigated for general neutralino candidates. The low 
B-ball decay temperature may also have interesting cosmological implications for the 
QCD phase transition and nucleosynthesis. Ultimately, we would hope to be able 
to provide a completely consistent model of B-ball Baryogenesis and the baryon to 
dark matter ratio in the context of a realistic SUSY infiation model which can satisfy 
the low reheating temperature constraint. Such a model would provide a basis for 
a complete cosmological scenario which would be able to naturally account for all 
cosmological observations (baryon number, dark matter and primordial fluctuations) 
and which would require no low energy physics beyond that of the MSSM. 
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Appendix. Some aspects of B-ball solutions 

In this Appendix we discuss some properties of numerical B-ball solutions and show 
that, for the purposes of discussing the thermalization and decay of the B-balls, a 
Gaussian Ansatz is a reasonable approximation to the exact B-ball solution. 

From the point of view of cosmology and phenomenology, the important quantites 
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are the energy and radius of the B-ball as a function of its charge B. The B-ball 
solution is of the form 



V2" 



(96) 



The energy and charge of the B-ball are then given by ||T0|, |T2 

2 



E 



(fx 
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and 



B 



(9^ 



where we rescale B such that B = +1 for the scalars. The equation of motion for a 
B-ball of a fixed value of u is given by 

2 , dU{(j)) 



r 



(99) 



where 0' = d(j)/dr. We require a solution such that 0(0) = ^'(0) = and as 
r 00. This corresponds to a tunnelling solution for the potential -~U{(f)) |jlO|, where 



(100) 



In practice, when obtaining numerical solutions, we vary 0(0) with these boundary 
conditions until the correct form of solution is obtained for a given u. The energy and 
charge of the solution are then calculated using the above expressions. 
For the d=6 direction, the B-ball equation is given by 
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where Uo is defined by 



uo^ - (1 + K) 



(102) 



and M is of the order of the scalar field within the B-ball. 

For the case of thin- walled B-balls, the initial value of is very close to 0ci, the value 
of for which the right-hand side of Eq. ( |101| ) vanishes. In this case will remain close 
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to 0ci, up to a radius of the order of 0J~Hog {(f)ci/S (f)(0)) where S(f){0) — (0ci — 0(0))- 
It will then decrease to zero over a distance Sr ~ o;"^, corresponding to the width of 
the wall of the B-ball. The radius of the thin- walled B-ball can be made arbitrarily 
large by choosing (50(0) small enough. 

For the case of thick-walled B-balls, the initial value of can be much smaller 
than 0ci- In this case the non-renormalizible terms may be neglected. In general, the 
right-hand side (RHS) of the B-ball equation vanishes for three values of (p, which 
correspond to (pci, <Pc2 and zero. 0c2 corresponds to the point at which, assuming that 
the non-renormalizible terms can be neglected, the first two terms on the RHS cancel, 

= (2^) ■ (^"^^ 

0c2 is an attractor, in the sense that if 0(0) is close to 0c2 it will tend towards 0c2 as r 
increases. We can now understand qualitatively how the thick-walled B-ball solution 
works. Suppose that initially 0(0) is large compared with 0c2 but small compared with 
(pel- Then as r begins to increase the RHS of the B-ball equation will be approximately 
—ojgCi)- Thus for small enough r the solution of the B-ball equation is 

Uo r 

This solution would become negative for r ^ u!~^. However, as approaches 0c2, the 
attractive nature of this value of alters the solution and, for one particular value of 
0(0) for a given ujo, causes the solution to "level out" and smoothly tend to zero as 
r — > 00. This gives the thick-wall B-ball solution for a given uJo- 

Numerically, we find the following properties for the thick-walled B-ball solutions; 

where Rc is defined as the radius within which 90% of the B-ball energy is found, and 

Uo ^ [3 - 4:)\K\^/^m . (106) 

Since typically is small compared with 1, we see that a; ~ m. 

We next show that a Gaussian Ansatz is a physically reasonable appoximation to 
the thick-wall B-ball solution. If we insert the Gaussian Ansatz, 

0(r) = 0(O)e":^ , (107) 
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in the left-hand side (LHS) of the B-ball equation, we obtain 



Inserting in the RHS gives 

- u;l<P - j r'cp ^{A + Br')<P , (109) 

where wc have set M = 0(0). Thus we see that the same form is obtained on the LHS 
and RHS. The exponential factor in the Gaussian will be correct up to a factor of the 
order of 1 if the values of A and B can be consistently approximated by a single value 
of R. (r as a function of cj) will be given correctly up to a factor of the order of the 
square root of the factor between the LHS and RHS). This requires that ~ 3|i^|m^ 
and ~ 2(|iir|m^)~^ in the numerical solution. 

Let us note some other properties of the Gaussian B-ball Ansatz. The total charge 
of the B-ball is given by 

3/2 



B 



/2r'^ /7r\ ' 

dr 47rr^u(t)le~'^ = (-) uj(j)lR^ . (110) 



The total energy of the B-ball is given by 

where the second term is the combined contribution from the potential energy and the 
charge term and we have used uj"^ ~ m^. Since R is large compared with for small 
\K\, the potential plus charge term dominates the energy. The radius within which 
90% of the energy is found, Rc, is then given by Rc — 1.2bR. The energy per unit 
charge is given by 

E / 3\K\\ 

where we have used the Gaussian result cu^ — 3|ir|m^. Although the energy per unit 
charge is larger than m, the mass of the scalar at small values of cf) will be of the form 
m(l + Q;|i^'|) (with a ^ 1) once the logarithmic correction to the potential is included, 
so that the binding energy per unit charge will be positive and of the order of l/^lm. 
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Comparing the Gaussian solution with the numerical solution, using the numerical 
value of Rc to normalize the Gaussian solution, we obtain R = kR{\K\^/'^m)-'^ with 
kn ~ 1.2, in good agreement with the value expected from the Gaussian Ansatz. 
The value of ujo from the numerical solution is about a factor of two larger than that 
expected from the Gaussian Ansatz. Thus we expect that r as a function of (p will be 
correctly given by the Gaussian Ansatz up to a factor of about 2. 

Physically, a factor of 2 variation in cUo correspond to a variation in or m by 

a factor of 2 in the real solution, which will not greatly alter the physical properties of 
the corresponding B-ball. Thus we expect the Gaussian to be a physically reasonable 
approximation to the B-ball solution. 
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